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Abstract 

We establish the incompressible Navier-Stokes limit for the discrete velocity model of 
the Boltzmann equation in any dimension of the physical space, for densities which re- 
main in a suitable small neighborhood of the global Maxwellian. Appropriately scaled 
families solutions of discrete Boltzmann equation are shown to have fluctuations that 
locally in time converge strongly to a limit governed by a solution of Incompressible 
Navier-Stokes provided that the initial fluctuation is smooth, and converges to appro- 
priate initial data. As applications of our results, we study the Carleman model and 
the one-dimensional Broadwcll model. 



1 Introduction and Main results 

> 

£T) . 1.1 Introduction 

o 

\Q ' The endeavor to understand how fluid dynamical equations can be derived from kinetic 

theory goes back to the founding works of Maxwell [38] and Boltzmann [16] . Most of deriva- 
tions are well understood at several formal levels by now, and yet their full mathematical 
justifications are still missing. Here we establish a so-called incompressible Navier-Stokes 
fluid dynamical limit for the discrete Boltzmann equation which is an evolution model for 
a gas which can attain only a finite number of velocities, one has to mention that Inoue 
and Nishida [30], and Caflish and Papanicolaou [19] have developed an asymptotic theory 
for the six velocity Broadwell's model. Further Studies have been developed for the two 
velocity Carleman's model by Kurtz [33] and by Lions and Toscani [36]. The interested 
reader can recover the pertinent literature in the survey by Lachowicz [34]. 

As known, see Bonilla and Soler [14], different asymptotic expansions can be developed 
for kinetic models thus obtaining various hydrodynamic descriptions. It has been recently 
recognized by several mathematicians, Bardos, Golse and Levermore [3, 4], Demasi, Es- 
posito and Lebowitz [21], that the incompressible Navier-Stokes equation can be obtained 
as the limit of the Boltzmann equation, when the Knudsen number e > go to zero. The 
validity of this fluid-dynamical approximation have been studied by several authors, see [2, 
5, 6, 10, 21, 25, 26, 37]. On the other hand, there is non proof concerning the convergence 
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of discrete velocity kinetic models to a solution of the incompressible Navier-Stokes equa- 
tions. However, although a complete proof is not given, some partial results may be found 
in [9, 10, 11, 28] in particular for Carleman and BroadwelPs models. In the present paper 
a rigorous proof of the connection between these two points of view is done, precising and 
completing some results which were announced in [4]. 

Let e > be the Knudsen number, and consider the scaled discrete Boltzmann equation 

ed t F i + v i -V x F i = -Q i {F,F), i = l,...,m, (1.1) 

where Fi = Fi(t, x) represents the mass density of gas particles with the constant velocity 
vector Vi = (vj, . . . ,vf) G W 1 at time t > and position x = (x±, . . . ,x n ) G W 1 ; Qi is 
a quadratic operator related to the binary collisions: 

Qi(F, G) = ^-J2 \ A k\( F kGi + FiGk) - Afj(FiGj + F 3 G t )\ , (1.2) 
lai jki 

where a>i are positive constants, and the terms Af- are the so-called transition rates refer- 
ring to the collisions 

(Vi,Vj) <-> (v k ,vi). 

The transition rates are positive constants which, according to the indistinguishability 
property of the gas particles and the reversibility of the collisions, fulfill the following 
properties: 

A*i — AV — A3i _ . . . _ \kl (-, o\ 

A lk — A kl — A kl — — l 1 -^ 

A detailed computation of the terms A\ l - can be performed by specializing the velocity 
discretization and analyzing the related collision mechanics. 

The above mathematical structure includes, as particular applications, various math- 
ematical models of discrete Boltzmann equation which can be recovered in the pertinent 
literature [18, 23]. It also refers to models with arbitrarily large number of velocities which 
have been studied in a number of recent papers, see among others [1, 27], considering their 
relatively simpler structure, with respect to the full Boltzmann equation, which may be 
exploited for scientific computing. Indeed, convergence to the full Boltzmann equation 
can be studied for models with arbitrary large number of velocities. 

Following [7, 18, 23], we shall introduce the basic concepts concerning (1.1) and sum- 
marize their properties which will be used later. 

Definition 1. A vector 4> = tu lim ..,0 m ) G M m is called a summational invariant if 

Kt \ai ctj a k ai J 

for all i, j, k, I = 1, . . . , m. 

We denote by M the set of summational invariants. Then < dimM < m because 
*(ai,...,a m ) £ M and M / W 71 . The following three conditions are equivalent: 
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(i) 4> £ M, 

(ii) ((f), Q(F, G)} = for all F,Ge R m , 

(iii) {<p, Q(F, F)) = for all F £ R m . 

Here (•, •) denotes the standard inner product in IR m . For the proof see [7, 18, 23]. 
Let F = t(Fi,...,F m ) G We write F > if F< > for all t = 1, . . . , m. 

Definition 2. A vector F = t(Fi,...,F m ) > is called a local Maxwellian if 
A%(FiFj - F k F t ) = 0, for all i,j,k,l = l,..., m. 

In particular, F > is called an absolute Maxwellian if it is a locally Maxwellian state 
and is independent of t and x. 

Consider the initial value problem for Eq. (1.1): 

n 

ed t F £ + V>d Xj F e = -Q(F £ ,F £ ), t > 0, xeR n , 

F £ (0,x)=F (x), (1.4) 

where 

V j = diag(uj,...,uf) , j = 1, . . . , n, 

and 

Q(F, G) = i(Q 1 (F,G),...,Q m (F,G))- 

According to [4], the discrete velocity Boltzmann equation (1.1) gives the incompressible 
Navier-Stokes equation in the limit e — > if F £ remains near an absolute Maxwellian M 
with the distance of order e. 

Let M > be an absolute Maxwellian state, and let 



We shall see the solution in the form 

F e (t,x) =M + eldf £ . 
Then Problem (1.4) can be transformed into the following one: 

1 n 11 

dtfe + - V Vid Xj f e + -Lf £ = -r(/ £ , f £ ), 

f £ (t = 0,x) = f (x) = A-t(F (x) - M), (1.5) 

where the operators L and T are given by 

Lf = -2A- 5 Q(M, A3 /) , (1.6) 
r(/, 5 ) = A-^Q(A§/,Al 5 ), (1.7) 

and have (see [7, 18, 23]) the following properties: 
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(i) L is real symmetric and positive semi-definite. It is null space is given by N(L) = 
A2M. We denote by {ej, i = 1, . . . , d} the orthonormal basis for N(L). 

(ii) T is bi-linear and satisfies T(f,g) G N(L) 1 - for any /, g E R m , where N(L) 1 - denotes 
the orthogonal complement of N(L) in M m . 

Existence of global solutions for discrete kinetic models have been obtained by several 
authors (see [15, 17] and the review article of Bellomo and Gustafsson [8]). However 
these results cannot be used here because they are based on the dispersion properties 
of the linearized equations and concern small (with respect to the collision operator) 
perturbations of the vacuum. It is easy to prove that any finite velocity model has a unique 
smooth solution during a finite time that depends on the size of the initial data and on 
the collision operator (or a global existence see [31, 32]). To the best of our knowledge it 
has not been proven that the time of existence of this smooth solution is independent of e. 

Now we state our main results. 



1.2 Main results 

Our main result is an existence theorem that holds for all e > and a proof of the 
validity of the fluid-dynamical approximation. To state our result precisely, one needs 
some function spaces. 

Let C(fl, X) and L°°(f2, X) denote the spaces of the continuous and bounded functions 
on 0, C W 1 with values in a Banach space X, respectively. 

H l denotes the I? (R n )-Sobolev space of order /, with the norm || • 

Theorem 1. Let (1.3) be assumed. Let n > 1 and I > § . If fo € H l (R n ), then there 
exists a positive constants To and k (depending only on \\fo\\i) such that the initial value 
problem (1.5) has a unique solution f £ G L°° ([0, T ],H l ) DC ([0, T^H 1 ' 1 (W 1 )) satisfying 



\\fe(t)\\l<k, (1.8) 

forte [0,T ]. 

If, in addition the initial data satisfies: 

/ e (0) = h £ + e 2 k e where k e £ H l and, 

h £ G N(L), \\T{h e , h e )\\i-! < Ce, H^Mz-i < Ce, 

lim \\h £ - hlU-! = 0. (1.9) 

one gets the strong convergence for the discrete Boltzmann equation to Incompressible 
Navier-Stokes equations: 

Theorem 2. Let f £ be as in Theorem 1. Then, as e — > 0, f £ — > / weakly * in 
L°° ([0, T],H l ) and strongly in C ([0, T], H 1 ' 1 ) for any T > 0, and the limit has the form 

d 

f = Y.p^ ( L1 °) 

i=i 
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where pi = (e^,/) satisfy: 

V x ■ (Ve <g> e • p) = 0. (1.11) 
Moreover, p(t, x) is a weak solution of the equations: 

f t - V 2 : {(L-\P^{Ve)) ® (Ve)) • p) + \v • (P^(y e )^=^) = (V ■ Ve (8) e) ■ 7r, 
p(t = 0) = (M, (1.12) 

2 / 2 2 \ 

where f = denotes the vector I , = . , . . . , , gn , . ) and £/te vectors 
V 2 :^®!-^? 1 ^)))/) and V • (P^Fe), 



VaM 



are <?roen fry: 




V 2 : (y e 0L- 1 (P ± (Fe)))p= 2E^^^' e ' 

\i=i fc=i 

Let (p £ )j = (ej, / e ). Since {e^, i = 1, . . . , d} forms an orthogonal system, p in (1.10) is 
given by pi = (ej, /). One gets 

Theorem 3. Xei &e assumed. Then, as e — > 0, p £ — > p weakly * in L°° ([0, T],ii~') 

and strongly in C ([0, T], id' -1 ) for any T > 0, and t/ie Zirmi /9 satisfies the incompressible 
Navier-Stokes equations (1.11)-(1.12). 

The system (1.11)— (1.12) are generalizations of the Navier-Stokes equations. The con- 
dition (1.11) generalizes the V • u = 0, V(p + 0) = conditions that arise in the classical 
incompressible Navier-Stokes limit (see [3]). tt is the lagrange multiplier corresponding to 
the constraint (1.11) and is the generalization of the classical pressure term. The quadratic 
term in (1.12) corresponds to the classical convection terms. The justification of these 
formal approximations for the classical Boltzmann equation has proven difficult because 
many regularity questions remain open for both these fluid systems as the Boltzmann 
equation. Two approaches to circumventing these difficulties have emerged recently. First 
some authors have studied direct derivations of linear or weakly nonlinear fluid dynamical 
systems, such as incompressible Navier-Stokes [21]. Their result requires smooth initial 
data and holds for as long as the limiting solution of the incompressible Navier-Stokes 
system is smooth. Second, some authors have abandoned the traditional expansion-based 
derivations in favor of moments based formal derivation [5, 6, 26]. In [6] it is shown 
that the solution of the Boltzmann equation considered over space of dimension three or 
more will be smooth for all time, with small initial data and will converge strongly to 
the solution of the incompressible Navier-Stokes equations. Recently the Incompressible 
Navier-Stokes fluid dynamical limit for the classical Boltzmann equation is considered 
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in [26] . It was shown that the scaled families of DiPerna-Lions renormalized solutions [22] 
have fluctuations that globally in time converge weakly to a limit governed by a solution 
of the incompressible Navier-Stokes equations due to Leray [35] provided that 

H(fe(0)) = 11 (MO) log / e (0) - / £ (0) + 1)M dvdx < Ce 2 , 

J Jl'xR 3 

and the fluid moments of their initial fluctuations converge to appropriate L 2 initial data. 
The proof use the averaging lemma (cf. [24]). The averaging lemma is valid for contin- 
uous solutions and has no counterpart for discrete velocity models (except in one space 
dimension (cf. Tartar [39])). However, in the proof of Theorem 2, one needs to evaluate 
the limit of the nonlinear moment 

(L- l {P^{Ve)\T{f e J e )). (1.13) 

Some uniform regularity estimates would likely be needed for obtaining the limit of non- 
linear terms (1.13). This term disappears in the case of Stokes limit (see [12]). 

Here we establish a so-called incompressible Navier-Stokes fluid dynamical limit (1.11)— 
(1.12) for the discrete Boltzmann equation in any dimension of the physical space. We 
solve (1.5) by using the principle of contraction mappings, we introduce the iteration 
scheme (Lemma 1 and Lemma 2) similar to one used in [19] to obtain the uniform estimate 
of the remainder term in the Hilbert expansion in order to justify the compressible Euler 
system from the one dimensionnal Broadwell. The convergence of the scheme to a solution 
of equation (1.5) is proved by using some properties of the operators L and T. The strong 
convergence of the solution of equation (1.5) as e — ► is proved by the uniform estimate 
and the equicontinuity in t € [0, T] of the solution with respect to £ £ (0, 1) (Lemma 3) 
provided that the initial fluctuation is smmoth, close to an N(L) element which converges 
to appropriate initial data. 

Remark 1. Put h £ = A2i( Ql ,a 2 ,...,a m )We, one gets 

T(h e ,h e ) = w 2 A~^Q(M,M) = 0. (1.14) 
An example of assumptions (1.9) can be given by: 

\\d x w £ \\i-\<Ce, lim \\w £ — w\\i-\ = 0. (1-15) 

£-+0 

h 2 

Remark 2. In the case where ^= G N(L), one has T(h £ ,h £ ) = (see Lemma 5). 
Therefore an other example can be given by: 

/ e (0) = h £ + e 2 k £ , where k £ € H l and, 
h 2 

h £ eN(L), ~~=== € N(L), ||aA||j_i < Ce, lim \\h £ - h]^ = 0. 

VaM e^o 

Remark 3. In [12] it was proven that the kinetic models (1.1) converge weakly and 
strongly to the linearized incompressible Navier-Stokes. The assumption T(h £ , h £ ) = 0(e) 
in H 1 ^ 1 is not necessary in this case and the time of existence of solutions goes to infinity 
as e — > 0. 
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Remark 4. The present work improves upon the result given in [12] for the Carleman and 
Broadwell's model, without assuming some restriction upon the scaling of the fluctuation 
with respect to Knudsen number, namely that the fluctuations should only be required to 
be of an order equal to the Knudsen number, but at the cost of being restricted to initial 
fluctuations (assuming T(h £ ,h £ ) = 0(e) in H 1-1 ) for the Broadwell model. 

The plan of this paper is as follows. Section 1 deals with the above introduction and 
main results. In Section 2, the formulation of the problem and the proof of uniform ex- 
istence Theorem are given. In Section 3, the estimate for dtf £ is shown. This estimate 
is used to prove the strong convergence of the solution to the solution of the non-linear 
incompressible Navier-Stokes equation. Finally, Sections 4 and 5 contains some appli- 
cations. As applications of our results, we are dealt with the Carleman model and the 
one-dimensional Broadwell model. 



2 Uniform existence 

To prove the existence of local solutions to (1.5), one has to get suitable "apriori" estimate. 
2.1 Estimates 

Lemma 1. Let z{t,x) be a given function of t and x such that, 

||z(t)||j<fc, (2.1) 
and let f(t,x) satisfy the linear system 

d t f + - £ J2v j d Xj f + ^Lf = - £ F(zJl 
i=i 

/(0,x) = / (x). (2.2) 
Then there exist To such that a constant k can be chosen such that 

sup 11/(^11, < k. (2.3) 

0<t<T 

Proof. From the theory of linear hyperbolic systems we know that (2.2) has a unique 
solution in L°° ([0,T],H l ) with f in L°° ([0,T],L 2 ). X-Mozilla-Status: 0000 
Taking the Fourier transform of (2.2) in x yields 

dtf + V j itjf + ^Lf = h(z, f). (2.4) 
i=i 

n 

Take the inner product (in C m ) of (2.4) with /. Since ^ V 3 Q and L are real symmetric, 

j i 

the real part of (2.4) is 

^f- + 1(L/, /) = \ Re (f(z, /), /), (2.5) 
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where (-, ■) denotes the standard inner product in C m . 

Noting that L is positive semi-definite and T 6 N(L)- 1 -, one obtains due to (2.5) and to 
the the inequality ab < \ (a 2 + fe 2 ) , the following estimate: 

^ + |l^/T < |f( ^" P±/l < ^rlf (z, /)| 2 + g^/l 2 , (2-6) 

where C\ is a constant and P 1 - is the orthogonal projection onto N(L) L . 
In particular (2.6) implies that 

^f<^r\t^f)\ 2 - (2-T) 

Noting that for Z > §, 

||r(/,0)||,<C 2 ||/||,||<7||i V (2.8) 

So by multiplying (2.7) by (l + |C| 2 )\ integrating over [0,t] x R" and using (2.8). One 
obtains from Plancherel's Theorem, the following inequality: 



< ll/ollf + ^ f ||/(*)||?||*(-)||, 2 da < ||/o|| 2 + Ck 2 T sup ||/(t)|| 2 , (2.9) 
b i jo te[o,T\ 



c 2 

with C = -jr-. 

Let Tq and A; are such that 



1 , l-Vl-4VCT||/o||, 



16C||/o|| 2 ' 2VCT 
for Te [0,T ]. 

Then the desired estimate (2.3) is an immediate consequence of (2.9). Thus the proof 
of Lemma 1 is completed. ■ 

We shall solve (1.5) by using Lemma 1 and the principle of contraction mappings. 
Define the iteration schemej/^} by 



fs — /0) 



N 



dtf" +1 + \ E V^f™ + ^/f +1 = iff, fs N+1 ) , 

f £ N+1 (0,x) = f (x), iV = 0,1,2... (2.10) 



Lemma 2. Let /o G H l . Then suitable constants Tq, k and A (A < 1) exist such that for 
any e > and for any t € [0,Tq], the following estimates are satisfied: 



fe N+1 \ 



|,<fc, (2.11) 
and 

f?l<C \l. (2.12) 
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Proof. Since \\f \\i = \\fo\\i < k, (2.11) follows thanks to Lemma 1. 
Let h£ = f£* +1 — /<f , therefore satisfy the equation: 

Applying the technique used in the proof of Lemma 1 one gets: 



N 



< 



1 

c 



r{h?-\f? +1 



+ 



Multiplying by (l + |C| 2 )\ integrating over [0, T] x RJ? and using (2.8). One can use again 
Plancherel's Theorem (2.11) to deduce that: 



c 2 

with C = -jr-. 

Since Cfc 2 T < 1, it follows 

,, N ..2 k 2 cT ||,JV-1||2 



(2.13) 



Put 



A 



A; 2 cT 



1-PcT' 

It is clear that A < 1 and (2.12) follows from (2.13). ■ 
2.2 Proof of Theorem 1 

In view of Lemma 2 the estimates (2.11), (2.12) imply that for each e > 0, {f^} is 
a Cauchy sequence in L°° ([0,T],H l ). Let denote its limit by f s (t) and note that it 
satisfies the estimate (1.8), i.e., this limit in L°° ([0, T], 

From (2.10) one sees that d t f^ +1 can be expressed in terms of sequences converging in 
L°° ([0, T^H 1 ' 1 ) as iV — ► +oo. The limit is 

1 n 11 

H £ = —J2 Vid Xj f £ - -Lf £ + -r(/ £ , f £ ). 

3=1 

Now let ty(t,x) be a C°° function of compact support in [0, T] x K. We have just seen 
that 

I {i{j(t,x),d t f £ N+1 )dtdx^ [ (il>(t,x),H e (t,x))dtdx, 

J[0,T]xK J[0,T]xK 

as N — > +oo. However 

/ Mt, x),d t f £ N+1 ) dtdx = - [ {dtW, x),f £ N+1 ) dtdx 

J[0,T]xK J[0,T]xK 

(d t ip(t,x),f £ )dtdx, 



[0,T]xK 



The Incompressible Navier-Stokes for the Non Linear Discrete Velocity Models 



435 



as N — > +00. Therefore H e (t) is identified with the distributions derivative in t of f £ . It 
follows that f £ satisfies the equation (1.5) and moreover ^ € L°° ([0, T], H 1 ^ 1 ); hence 
f e G C([0,r],fl' z - 1 ). X-Mozilla-Status: 0000 



3 Hydrodynamical limit 

This section deals with the strong convergence of the solution of the discrete velocity 
model (1.1) towards solution of the nonlinear incompressible Navier-Stokes equations 
(1.11)— (1.12). In order to derive the hydrodynamical limit, some uniform regularity esti- 
mates would likely be needed for obtaining the limit of the nonlinear terms. 



3.1 Strong convergence 

To prove Theorem 2, one needs more than the uniform bound (1.8) for f £ . For this purpose 
we assume that hypothesis (1.9) holds and we show the uniform bound for 
The uniform equicontinuity in t is given by the following: 

Lemma 3. Assume that hypothesis (1.9) holds and let I > ^ + 1. Then 

dfe 



( Ck 2 T\ 

<Cexp(__J, Vte[0,T], e€(0,l), 



dt 

where the constant C does not depend on e. 



(3.1) 



Proof. By differentiating the equation (1.5) with respect to t and by taking the Fourier 
transform of the equation obtained, one has 

j t dtfe + -Y^i^dtU + j 2 Ld t f e = - £ d t t(f £ J £ ). (3.2) 

Taking the inner product (in C m ) of (3.2) with d t f £ . The real part of the resulting equality 
is 

+ \(LdJ £ ,dJ) = -Re{f(d t f,f),dtf.) 
Since L is positive semi-definite and V e N(L)- 1 -, we get 

dt\dtfe\ 2 , C 1|D± |f(d t / £ ,/ £ )||P^/ £ | 



+ -±\P^ dt f £ \< < 



1 . .s. 2 Ci ± f 2 



So, multiplying by (l + |C| 2 )' \ integrating over and using (2.8) one gets 

W/elli 2 -i < ^l|r(ft/e,/ e )ll, 2 _i < c\\d t f £ \\U\\f £ \\U. (3.3) 
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Using Gronwall's inequality, one reduces (3.3) to 

'Ck 2 T x 



\d t fs\\i-i < exp 



\dtfe(0)\\ 



i-i- 



By using equation (1.5) in order to express dtf £ (0) in terms of the initial data we have 

'Ck 2 T" 



\d t fe\\l-i < exp 



3=1 



+ ^||L/ ||,_i + i||r(/o,/o)||i-i 



l-l 



Taking into account hypothesis (1.9), one gets 

'Ck 2 T\ (\ 



%fe\\i-i < Cexp 



d x h E \\i-i + £\\d x k £ \\i-i + ||Lfe e |b_i 



1 



+ £\\h £ \\i\\k £ \\i + £ d \\k e \\f + -||r(h e ,/i e )||i_i ) < Cexp 



Ck 2 T 



The proof of (3.1) is now complete. ■ 

Conclusion. From Lemma 3 we conclude the following: the solution f £ is bounded in 
C([0,T], ii"'- 1 ), uniformly for e > 0, for t in any compact subset of the interval [0, T], 
moreover f £ satisfies the bound (3.1). Therefore by the Ascoli-Arzela lemma we can 
choose a convergent subsequence f £j , where £j — > such that 



/e, > f 



m 



c([0,T\,H l - i y 



(3.4) 



and the limit function satisfies the bound (1.6). 

Thanks to convergence results (3.4), the estimate (1.8) allow us to deduce (possibly 
taking subsequences) the following convergence 

r(/ £ .,/ £ .)^r(/,/) strongly in H 1 ' 1 . 



(3.5) 



3.2 Passage to the limit 



This subsection is devoted to the connection between the discrete velocity kinetic equa- 
tions (1.1) and the incompressible Navier-Stokes equations (1.11)— (1.12), we use the tech- 
niques described in [3]. In order to get these equations, one needs the following: 



Lemma 4. Let f,g € M m . One gets 

(/, Lg) = \ E ^(M t M + M k M{)Ut + /* - ft ~ ftM + 9* - 9* k - 9*i 



(3.6) 



i,j,k,l 



where 



fi 



Proof. See [7] and [23]. 



i = 1, . . . , m. 
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Lemma 5. Let f £ N(L). One has 

r</ ' /)= Hiy- 

Proof. Using (1.7), one gets 



(3.7) 



j,k,l 



ST I & ^ M k M l f f A kl v , f 



a.;* a. 



(3.8) 



On the other hand, if / € N(L), so one gets by (3.6) the following identity: 

fi fk 



A ki 



fi _|_ fj 



VotiMi ^oijMj s/cqMi Va k M k/ 
Thus, we have 



0, Vi,j,k,l. 



2^" 



fifj 



A kl 



f 2 f 2 
h _|_ Ik 



ft ft 



^/aiajMiMj 
Hence substituting this identity into (3.8) yields: 



+ 



2flfk 



a t Mi a k M k a { Mi ajMj ^aia k M t Mk 



r(/,/)= J>2 

V j,k,l 



-fkfl - 



--fkfl ~ 7,—^= TTJl 



^aiaia k Mi yJa,iaia k M k Mi 2 JaictiMi ' 



1 y/WiMj f2 t 1 ^MlMj 2 1 ^Wj 2 

Ji ' o / — J 



2 yJa~iOi k Mk J K ' 2 aiy/oiMi Jt ' Ictj^fcTi 3 



-fk + o" 



(3.9) 



Noting that 

A%(MiMj - M k Mt) = 0, V i, j, k, I. 

So, 



A kl 



f k fi 



^a i aia k M l 
which with (3.9) imply 

r(/,/) = ~|E45 



- fkfl = 0, 



y/aiaia k M k Mi' 



\j,k,l 



fttiMj 2 y/ttiMj 2 
-J) 



faiaiMi 



ft + 



— ft 

a~ia k M k JK a>iyfa~iMi Jl aj^/ai 3 



M fc ,/f , .U; /I Mj / 2 M, ft 



+ 



yJouMi an yJaiMi a k y/oLiMi oii \JctiMi atj 
Therefore the desired identity (3.7) follows from (3.10) and (1.6). 



. (3.10) 
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Lemma 6. Let x be a test function such that 

n d 

J2J2 d *iXr(er,Vi ei ) = 0, l = l,...,d. (3.11) 

j=l r =l 

Then 

n 

(J2dx j (fe,P(V j (e l )), X ) = 0, l = l,...,d, (3.12) 

where P is the projection onto N(L). 
Proof. One has 

n n 

(Y,d Xj (f e ,P(Vi( ei )), X ) = (Z&APf'^W)'*) 

3=1 3=1 

n d n d d 

= ^2( d ^j(^2 a e,re r ,V J (ei)),x) = ^2^2^2(d Xj a £:r (e r ,V 3 {ei)},xi}- 

j=l r=l j=l r=l 1=1 

Integrating by parts and using (3.5), one gets 

n n d d 

£>.</ £ ,P(W(eO),X> = -EEE<a e ,r<er,W(e,)>,a Xj . X i> = 0. 

3=1 j=l r=l 1=1 

This ends the proof of Lemma 6. ■ 

Completion of the proof. Multiplying the equation (1.5) by e 2 , letting e go to zero 
and using (3.4), yields the relation 

Lf = 0. 

This implies that / G N(L) and thus can be written according to the formula (1.10). 
The scaled local conservation laws associated with the kinetic equation (1.5) are 

1 n 

d t (f E ,e i ) + -J29x j (V J f E ,e i )=0, i = l,...,d. (3.13) 
£ j=i 

Letting e tend to zero and applying formula (1.8) yields condition (1.11). 

As was done in [3] for the classical continuous Boltzmann equation, we divide by e and 
we decompose the flux in the form 

1 n 1 n 

d t (f £ ,e t ) + -^^,(/ £ ,P ± (^e i )) + -Y,d Xj {f e ,P{Viei)) = 0. (3.14) 

£ 3=1 £ 3=1 

Using 3.11 (Lemma 6), the second term of (3.14) can be eliminated upon integra- 
ting (3.14) against test functions satisfying (3.11). 
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The limit of the first flux term in (3.14) is computed using the fact that P^iV^ei) £ 
Ran(L) along with the kinetic equation (1.5) to obtain 

n n 

d t (f £ , ei ) + E d Xj {-edtf e - E vkd *Je + r(/ £ , f e ), L- 1 P ± (V^e i )) 
3=1 k=l 

1 n 

+ -^9 Xj {f e ,P(V^i)) = 0, i = l,...,d. (3.15) 

£ 3=1 

Letting e goes to zero, 

ran ra ra 

EE^a^v^-^^'^-e 

3=1 k=l j=l k=l 

= V 2 : (Ve®L- 1 (P ± (Ve)))p, (3.16) 

The limiting quadratic term of (3.15) must be evaluated further in order to bring it 
into the form that appears in (1.12). Since / € N(L) the convergence (3.5) and the 
identity (3.7) can be employed to show the following convergence in D' t x as e — > 0, 

ra n 

E d *i We, fe^L^P 1 - (We,)) - E 9 *i W, /), 

= if^M/^)) = I V • {P ± (Ve), ^=£). (3.17) 

Passing to the limit e — > in equation (3.15), the above convergence (3.16) and (3.17) 
suffices to write the limit equation (1.12). 

By integrating the equation (3.15) over t, we have 



{fe,ei)-{f e (0),ei) 



ft f n n 

= ~ / \^ dx ^- £d ^-T, Vkd ^ + nfeJe),L' 1 P ± (V^e l )} 
J ° \j=l k=l 
ra \ 

+ -Y J d Xj {f e ,P{V^e i )){s)\ ds = 0, i = l,...,d. 

Let e — > and putting t = 0, it follows 

</,e i )(0) = (h,e i ). (3.18) 



Therefore the desired initial condition for the system (1.12) follows from (3.18). 
In the next, we will apply the results discussed in Section 1 to two examples. 
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4 Example, I (Carlemann model) 

The simplest one-dimensional discrete- velocity models of the Botzmann equation are cer- 
tainly those with two velocities. These models describe the evolution of the velocity 
distribution of a fictitious gas composed of two kinds of particles that move parallel to 
the x-axis with constant speed equal to one, either in the positive x-direction with a den- 
sity F\, or in the negative x-direction with a density F2. The corresponding dimensionless 
hyperbolic system is given by (see [20]): 



dF? dF? 1 



e- 



1 , U1 1 _ x ( W e2 p£ 2\ 



dt dx e 
dFZ dFZ 1 



e 



2 U1 2 _ x / 77e2 



- (Ff 2 - F% 2 ) , 



dt dx 

Ff(t = 0,x)=F°, i = 1.2. (4.1) 
Here we have (ai, 02) = (1, 1) and 

^11 = ^22 = 1) and ^tz = otherwise. 

The condition (1.3) is satisfied. The space M of summational invariants consists of vectors 
<j) = £(0 1; 2 ) satisfying (f>\ — $2 = 0. Therefore M is spanned by (j) = (1, 1). 

On the other hand a locally Maxwellian state is a vector F = trp lt p 2 \ > satisfying 
Ff = Ff. Let M = (1,1) be an absolute Maxwellian state: M = (1,1). Set F(t,x) = 
M + £Aa/(t, x), here A = and substitute it into (4.1), we get: 

dtfe + \d x Vf £ + ^Lf £ = hr(f e , f e ), 

f £ (t = 0,x) =/ (x), (4.2) 



1 -n F a 



where 



-11/' V -1 

and 

r(/,/) = (/i 2 -/ 2 2 ) 

As e — > 0, we get that / £ — > / in (distribution sense) with f = p-e and p is solution 
of the heat equation: 

d t p = \dlp. (4.3) 

It follows from the uniqueness of the solution of the initial value problem for (4.3) that 
all sequences of f e as e — > give the same system (4.3) in the limit. 
Thus we have proved 

Theorem 4. Let /o £ i? 1 (IR), i/ien i/iere exists a positive constants Tq and k (de- 
pending only on \\fo\\i) such that the initial value problem (4-2) has a unique solution 
f £ e L°° ([0,T ],H l (R))nC{[0,T],L 2 (R)) satisfying 

\\fe(t)\\l < k, (4.4) 

forte [0,T ]. 
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If in addition /q satisfies 

/e(0) =t (heM +e 2 k £ , \\d x h £ \\ 2 < Ce, lim\\h £ - h\\ 2 = 0, (4.5) 
we get 

Theorem 5. As e — > 0, f e — > f = p ■ e weakly in L°° ([0, To],!! 1 ) and strongly in 
C ([0, T], L 2 (x)) and the limit satisfies the heat equation: 

dt P =\d 2 x p, p(0,x) = (h,e). (4.6) 

5 Example, II (Broadwell model) 

We shall first describe the six-velocity model gas considered by Broadwell [17] and in [29] 
and then specialize it to one dimension. We consider the following equation: 

+ u i ■ VF * = ^-(*i+i*i-H + F i+2 F i+5 - 2FiF i+3 ), i = 0,..., 5 (5.1) 
at 6e 

Here Fi are the densities of the particles with velocity 

Ui = ^ccos (6 + i^j , csin (j) + > 

and B denotes the collision frequency. 

We specialize (5.1) to solutions that do not depend on y, and 9 = 0, and such that 
F 5 = Fi and F A = F 2 . We can write (5.1) for F 1 , F 2 , F 3 and F 4 : 

dF dF _2B 

+ ~dx~~ 3F ( 1 2 " 3) ' 

Lit \JjJu UG 

F x IdFx B . „ „ . 

dF 2 ldF 2 B 

£ ^f ~ 2~dx~ ~ Te {FoF3 ~ FlF2) ' 
dF 3 dF 3 _2B 

£ ^f-^-^ {FlF2 - F ° F3) - (5 - 2) 

By the form of second member we have («o, u\,a 2 , a 3 ) = (^, 1, 1, |) and the relation 
(1,3) is satisfied. To get the equation (1.11)— (1.12) for the Broadwell's model we need some 
preparations. The space M of summational invariants consists of vectors <f> = i^,^,^,^) 
satisfying 2<j>\ + 204 = 02 + 03- Therefore M is spanned by {0i, cf> 2 , <j) 3 }, where 

01 = <(i,o,o,-i)' ^ = *(o,i,oi)> 03 = *(o,o,i,i) 

On the other hand a locally Maxwellian state is a vector F = t(F ,F 1 ,F 2 ,F a ) > satisfying 
F 2 Fi = FqF 3 . Let M = ini^i) be an absolute Maxwellian state. Set F(t,x) = M + 

hhf(t,x) for A = diag (2, 1, 1, 2) and substitute it into (5.2): 

ft + ^ + ^ = -T(f,f), 

e e z e 

f(t = 0,x)=f (x), (5.3) 
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where V = diag (l, 5, —5, — l) and 



3 



/ -2 ^2 \/2 -2 \ 

v/2 -1 -1 \/2 

v/2 -1 -1 V2 

\ -2 v/2 -2 / 

5 



r(/,/) - — (/2/1 - 2/0/3)*^ -i-i lV ^)- 



(5.4) 



(5.5) 



Since iV(L) = A2M, a simple calculations gives the orthonormal basis {ei,i = 1,2,3} for 
N(L): 



ei 



^,0,0,-^ 



e2 = 



_L iL — !— 



e3 



By letting e — > in (5.3) and using the results of Section 1, we get f e 
(distribution sense) with / = p • e, where 

P ~ t(pi,p2,P3)> e — Hei,e2,e3)) 



and from (1.9) we get 



P3 



2v/3 



P2 



Lemma 7. Let f = p • e satisfying (5.6), then 

r(/,/) = o, 

(Vf,h)=0 for any h € N(L). 



(5.6) 



(5.7) 
(5.8) 



Proof. By using (5.5) and (5.6), we get (5.7). (5.8) is obtained by a direct calculations. 
This completes the proof. ■ 



Lemma 8. Let f = p • e satisfying (5.6), then 



(L-\P^(V ei )),Vf} = ^, 
{L-\P ± (Ve 2 )),Vf) = 



Pi 



(5.9) 
(5.10) 

(5.11) 



Proof. The proof of this lemma follows from Lemma 7. ■ 

Using the results of Section 1 combined with Lemma 8, one gets: 

Theorem 6. Let fo £ if 1 (R), then there exists a positive constants To and k (de- 
pending only on \\fo\\i) such that the initial value problem (5.3) has a unique solution 
f e G L°° {[0,T },H l (R))nC{[0,T],L 2 (R)) satisfying 



ll/ e (*)lli<*> 
forte [0,T ]. 



(5.12) 
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If in addition the initial condition satisfies: 

/ e (0) = h e + e 2 k £ , Lh e = 0, HIX^A)^ < Ce, 

\\d x h £ \\ 2 < Ce, ]hn\\h e -h\\ 2 = 0, (5.13) 
one gets: 

Theorem 7. As e — > 0, f £ —> f = p • e weakly in L°° ([0, To], .ff 1 ) and strongly in 
C ([0, T], L 2 (x)) and the limit satisfies the heat equation: 

PlO " 7?P20 + 2 ^|p30 

Pi(* = 0) = — , (5.14) 

where 

p i0 = (h,ei), i = 1,2,3. 

Remark 5. One can remove the assumption ||r(/i e , h £ )\\2 < Ce in (5.13). A price to be 
paid for this improvement is that: we see the solution of (5.2) in the form F £ = M+e^>(e)/ £ 
with <P{e) = 0(e) (see [12]). 
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